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ABSTRACT 

We present a new analysis of the infrared (IR) emission from the ejecta of SN1987A covering days 
615, 775, 1144, 8515, and 9090 after the explosion. We show that the observations are consistent 
with the rapid formation of about 0.4 M 0 of dust, consisting of mostly silicates (MgSiCU), near day 
615, and evolving to about 0.45 Mq of composite dust grains consisting of ~ 0.4 Mq of silicates and 
~ 0.05 Mq of amorphous carbon after day ~ 8500. The proposed scenario challenges previous claims 
that dust in SN ejecta is predominantly carbon, and that it grew from an initial mass of ~ 10“ 3 Mq, to 
over 0.5 Mq by cold accretion. It alleviates several problems with previous interpretations of the data: 
(1) it reconciles the abundances of silicon, magnesium, and carbon with the upper limits imposed by 
nucleosynthesis calculations; (2) it eliminates the requirement that most of the dust observed around 
day 9000 grew by cold accretion onto the ~ 10 -3 M 0 of dust previously inferred for days 615 and 775 
after the explosion; and (3) establishes the dominance of silicate over carbon dust in the SN ejecta. 
At early epochs, the IR luminosity of the dust is powered by the radioactive decay of 56 Co, and at 
late times by at least (1.3 — 1.6) x 10 -4 M 0 of 44 Ti. Even if only a fraction > 10% of the silicate dust 
survives the injection into the ISM, the observations firmly establish the role of core collapse SNe as 
the major source of thermally condensed silicate dust in the universe. 

Subject headings: dust, extinction infrared: general - Magellanic Clouds — supernovae: individual 
(SN1987A) 


1. INTRODUCTION 

First evidence for dust emission from the ejecta of 
SN1987A was provided by high-resolution rnid-IR im¬ 
ages obtained with the Thermal Region Camera and 
Spectr ograp h (T -ReCS) on the Gemini South 8m tele¬ 
scope dBouchet et alj|2004lh Spatially resolved images 
of SN1987A at 450 and 870 /im obtained with the Ata¬ 
cama Large Milimmeter/Submillimeter Array (ALMA) 
between days ~ 9200 and ~ 9300 after the explo¬ 
sion provided further proof o f thermal IR emission 
from the supernova (SN) ejecta (llndebetouw et al1l20l4 
IZanardo et alJl2014jh Combined with almost contempo¬ 
raneous data obtained on day 9090 after the explosion 
with the Herschel satellite, the infrared (IR) observa¬ 
tions have been interpreted as evi dence for t he presence 
of a large amount of cold dust (|Matsuura et all 1201 lL 
120151) . If the emission arises from only silicates in the 
form of Mg 2 SiC >4 then the observations require 2.4 M 0 
of dust, that is, the precipitation of, respectively, 0.48 
and 0.8 Mq, of Si and Mg from the gas phase of the 
ejecta. This mass greatly exceed the SN yield of either 
Si or Mg (~ 0.1 M 0 , see Tableland references therein). 
The mass of Si required to fit the IR spectrum can be 
reduced by the inclusion of carbon dust. A combina¬ 
tion of silicate and amorphous carbon reduces the mass 
of silicates to 0.5 Mq. However, it requires an excessive 
amount of carbon, about (0.3-0.5) Mq or over 2 to 3 
times the yield of C, to be locked up in dust. 

The mass of dust inferred from the late time (t > 
8500 d) IR emission also greatly exceeds the ~ ICR 3 M 0 
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of dust derived from early obse rvations of SN1987 A 
on days 615, 775, and 1144 (iMoselev et all 119891 
iWooden et all 119931 : iDwek etldl 1199211 . This differ¬ 
ence has been interpreted as evidence for the contin¬ 
uous growth by accretion of refractory elements onto 
the seed nuclei that were initially formed in the ejecta 
(1 Wesson et all 120151 . and references therein). The ejecta 
temperat ure after day ~ 1200 h as dropped below 

~ 500 K dSarangi fc CherchnefHl2013l) . An accretion sce¬ 
nario for the evolution of the dust will therefore require 
most of the growth to take place by cold accretion, so 
that the final structure of the dust grains will primarily 
consist of a mantle, loosely bound with surface energies 
< 1 eV, surrounding a small thermally-condensed core. 
Such dust grains will not survive t he reverse shock ex - 
pected to propagate into the ejecta (jNozawa et al.ll2007lh 
contrary to the observed ~ 0.04 M 0 of dust that was sub¬ 
jected to and survived the reverse s hock i n Cas A or in 
the ~10,000 old Sgr East remnant (lArendt et all 12014 
IBarlow et all 1201(1 : Rho et akll2008t lLau et all 12015lb 

Also, the absence of the 9.7 and 18 fin i emission 
features in the IR spectra on days 615 and 775 has 
been interpreted as evidence that the newly -formed 
dust c onsisted of mo s tly car bon dust (IWooden et al.1 
119931 lErcolano et al.l I2007al lbt IWesson et al.1 120 1 5l) . 
Using the observed spectra of days 615 and 775 and 
a fully 3D radiative transfer model, IWesson et al.l 
(|2015l) concluded that about (4 — 8) x ICR 4 M 0 of 
carbon condensed in the ejecta of SN1987A by those 
epochs. In particular, they derived an upper limit of 
~ ICR 4 Mq on the mass of silicates dust that could 
have formed in the ejecta. This low yield is puzzling, 
considering the preponderance of Si and Mg in the 
ejecta, and that SNe should be the major source of 
thermally-condensed dust in the local and early universe 
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In this paper we present an alternative model for the 
evolution of dust in the ejecta of SN1987A, and show 
that the currently observed dust mass in the ejecta on 
days ~ 8500 — 9100 could already have been present at 
the ear lier epochs. The conclu sions of lErcolano et all 
(|2007af l and I Wesson et alJ (|2015t f represent just one pos¬ 
sible, albeit troublesome, scenario for the evolution of 
the dust mass in SN1987A. Large amounts of dust could 
have been hidden in opti cally thick regions of th e ejecta, a 
possib ility first raised bv lLucv et al.l (|1989f) and lDanzigerl 
(|2009f h and now supported by the ALMA observations. 
We emphasize that our scenario for the evolution of dust 
mass in SN1987A is not unique, but a viable alternative 
which addresses the difficulties in previous models: (1) 
the large mass of either Si, Mg, or C that needs to be 
locked up in dust to account for the ALMA observations; 
(2) the growth of dust mass by cold accretion, creating 
loosely bound mantles that will not survive the passage 
of the reverse shock; and (3) the absence of any signifi¬ 
cant amount of silicate dust in the ejecta. 

Our alternative model has important implications for 
the composition of dust in the early universe (z > 7), be¬ 
fore carbon-rich AGB stars left the main sequence. Even 
if only a fraction > 10% of the silicate dust survives 
the injection into the ISM, the observations of SN1987A 
firmly establish the role of core collapse SNe as the major 
source of thermally condensed silicate dust in the early 
universe. Interstellar dust is then silicate-rich, with im¬ 
portant observational consequences for galactic extinc- 
tion co rr ections an d photometric redshift determinations 
( Dwek et al.ll20lJ h 

The paper is organized as follows. We first present 
the basic equations for the escape of IR photons from a 
homogeneous dusty sphere, the dust model used in our 
fitting procedure (Section 2). In Section 3 we present the 
fits to the observations, and the evolution of dust mass. 
We first start with the analysis of the late-time emission 
at epochs 8515 and 9090, when the ejecta was optically 
thin, and construct a model in which the IR emission 
arises from elongated composite grains consisting of a 
silicate MgSiC >3 matrix with amorphous carbon (AC) in¬ 
clusions. Such configuration can arise as a results of the 
coagulation of silicates and AC grains at late times in the 
ejecta. The mass of Mg, Si, and C needed to be locked 
up in the dust is entirely consistent with nucleosynthesis 
constraints. We then show that the same mass of refrac¬ 
tory elements could have already condensed out of the 
gas at the early epochs (days 615, 775, and 1144). Here 
the calculations assumed that the ejecta is homogeneous. 
A clumpy ejecta model is described in Section 4. In Sec¬ 
tion 5 we calculate the evolution of the IR luminosity, 
and compare it to the energy input by radioactivity. A 
summary and the astrophysical implications of our re¬ 
sults are presented in Section 6. 

2. THE ESCAPE OF INFRARED PHOTONS FROM A 
HOMOGENEOUS DUSTY SPHERE 

The general equation for the flux emerging from a ho¬ 
mogeneous dusty sphere of radius R at distance D is 


given by: 

Su( A) = 4A ^ A) T d ) P esc [r{ A)] , (1) 

where M d is the mass of dust in the sphere, k(A) is the 
mass absorption coefficient of the dust, and B„{ A, Td) 
is the Planck function at wavelength A and dust tem¬ 
perature Td- The function P esc is the probability for 
an IR photon to escap e from the sphere and given by 
ICox fc Mathewsl (|1969f) : IQstcrbrock fc Ferlandl (j2006f l: 
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where r is the radial optical depth of the sphere given 
by: 
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when r << 1 (4) 
when r —> oo . 


The escape formalism explicitly takes the distribution 
of optical depth in the sphere into account. As r becomes 
very large, P esc (X) ~ 1/k(A) [eq. (0 and (0], with the 
consequence that the spectral signatures of the dust are 
cancelled in eq. ©• Furthermore, all information on the 
dust mass in the sphere becomes indeterminate. 

Equation m and 0 can be readily generalized to a 
dust population characterized by different compositions, 
masses, and temperatures. The flux is then given by the 
sum over dust components: 


s*( A) =£ 4M 4 ^ 2 (A) ^( A > T ^) Pesc[r( A)] (5) 
j 


where r is the optical depth in the presence of multiple 
dust components: 

? ( A ) = db Md ’j ^ ( A ) ( 0 ) 

j 


where M d j, Kj( A), and T d j are the masses, mass absorp¬ 
tion coefficients, and temperatures corresponding to the 
different grain compositions j. 

In the high optical depth limit, S v ( A) will still be a 
blackbody, as long as all dust components are at the 
same temperature. If the temperatures differ, then the 
emission features, generated by k(A), and the absorption 
features, generated by P esc ( A), will not cancel perfectly, 
and the total spectrum will still retain an imprint of the 
individual dust spectra. The former case is illustrated 
by our fit at Day 775, and the latter case by fits at Day 
615 and 1144, as shown later in Figures [3] and Q] and 
Tabled 


3. THE EVOLUTION OF DUST MASS 
3.1. Nucleosynthetic Constraints on the Dust Mass 

The mass of dust that can form in the SN ejecta is lim¬ 
ited by the mass of refractory elements produced during 
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Figure 1. The mass absorption coefficients used in this paper. Left panel: Spherical silicate (MgSiOs) and amorphous carbon grains; 
Right panel: Ellipsoidal composite grains. The blue curve corresponds to pure ellipsoidal MgSiOs grains. Successive curves represent 
the effect of added AC inclusions with fractional volume filling factor of 0.016, 0.028, 0.05, 0.09, 0.16, 0.28, and 0.50. The green curve 
represents the k used to fit the spectra of epochs 8515 and 9090, and corresponds to grains with MgSiC >3 and AC volume filling factors of, 
respectively, 82 and 18%. 


Table 1 

Elemental Yields of a 20 Mq Star 


Metallicity 

C 

Mg 

Si 

Fe 

Reference 1 

0.004 

0.020 

0.008 2 

0.097 

0.245 

0.144 

0.101 

0.095 

0.098 

0.124 

0.068 

0.096 

0.074 

0.089 

0.080 

[i] 

0.008 

0.250 

0.096 

0.11 

0.075 

[2] 


1 Refe rences: fll-N omoto et al.l 1120061 ): [2]- IWooslev et all 
ilT9Mf 

2 Yields for this metallicity were obtained by logarithmic 
interpolation between Z=0.004 and 0.02 


(ALMA) at far-IR wavelengths, and found to have 
an elliptical shape with dimensions corresponding to 
freely expansion velocities of 1350 ± 150 and 750 ± 
250 km s" 1 , along the ma jor and minor axes of the ellipse 
(llndebetouw et al.ll2014D . For simplicity we will adopt a 
spherical shape for the ejecta with a constant expansion 
velocity of v = (1350 x 750 2 ) 1 / 3 = 910 km s _1 . The 
radius of the ejecta at the different epochs is given in 
Table [2] On day 9000 it is about 7 x 10 16 cm, and its 
radial 200 /iin optical depth is ~ 0.3. Since the ejecta is 
optically thin at these epochs, we start our analysis with 
the spectra of those epochs. 


the quiescent and explosive phases of the evolution of the 
progenitor star. SN1987A is believed to be the end stage 
of the evolution of a ~ 20 M^ star, identified a s the 
blue supergiant Sk -69 202 (IGilmozzi et al.l lT987l j . Its 
initial metallicity is presumably that of the LMC, which 
is ~ 40% solar or Z=0.008. Table Q] presents the post ex¬ 
plosion elemental yields of the major refractory ele ments 
in the ejecta. The yields from lNomoto et ahl (|2006f) were 
interpolated between t he tabulated values f or Z=0.004 
and 0.02. The yields of lWooslev et al.l (I1988T ) were taken 
from his initial model with LMC abundances, that was 
calculated to fit the early SN light curve. Both mod¬ 
els predict that Mg and Si are produced in about equal 
amounts of ~ 0.1 Mg, but differ in the mass of predicted 
carbon, which varies between 0.144 and 0.25 Mg. Both 
models produce the right amount of Fe (synthesized as 
56 Ni) to power the SN light curve. The iron resides in 
a hot low density radioactive bubble, and is not likely 
to condense into dust at early times. We therefore only 
considered silicate and carbon dust in our models. Since 
the masses of Mg and Si are about equal, the simplest 
way to incorporate them in the dust is to assume that 
the dust consists of SiC >2 and MgO clusters or MgSiOs 
grains. In our models we assume that the silicates are 
in the form of MgSiOs, and the carbon in the form of 
amorphous carbon (AC) dust. 

3.2. Abundances Constrained Fits to the Infrared 
Spectra of SN1987A 

The inner ejecta of the SN was well resolved with 
the Atacama Large Millimeter/Submillimeter Array 


3.2.1. Days 8515 and 9090 

High resolution ALMA images at 450 and 850 /im 
clearly resolved the emissi on from the ejec ta and 
the circumstellar ring (llndebetouw et all 12014 
iZanardo et ahl I2014H. Sp e ctral observations, sum- 
marized bv iMatsuura et ahl (|2015lh showed that the 
far-IR/submm fluxes obtained by Spitzer and ALMA 
observations are dominated by continuum emission from 
ejecta dust, except for the 70 pm flux, which may be 
contaminated by the [O I] 63pm line and by the contin- 
uum emission from the hot dust in the circumstellar ring 
dDwek et al.l l2010h . The combined Spitzer and ALMA 
measurements of the far-IR and submillimeter flux 
densit ies of SN1987A used in our analysis were taken 
from (jMatsuura et al.l [2015L their Table 1). We adopt 
here their estimate that at least 2/3 of the observed flux 
in the 70 pm band arises from the ejecta dust. 

To minimize the mass of dust in the ejecta, 
IMatsuura et ah ( 2015 11 considered the possibility that the 
dust grains are in the shape of elongated ellipsoids. Elon¬ 
gated grains radiate more efficiently at long wavelengths, 
thereby lowering the mass needed to produce the emis¬ 
sion. However, the shape effect in AC grains does not 
only increase their long wavelength emissivity, it also 
considerably flattens it with wavelength. Since most of 
the emission in the Matsuura model is attributed to AC 
dust, the resulting spectrum became too broad to fit the 
observations. To circumvent these problems we adopt 
a dust model consisting of ellipsoidal composite grains, 
comprising mostly of a silicate (MgSiOs) matrix with 
AC inclusions. The elongated shape increased the dust 
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emissivity, and the addition of AC inclusions further in¬ 
creased the dust emissivity, with only minor flattening of 
its wavelength dependence. 

Figure [T] depicts the various mass absorption coeffi¬ 
cients considered in our models. In the Rayleigh limit, 
where a << A, the mass absorption coefficient of the 
dust is independent of grain radius. The left panel shows 
the wavelength dependence of k(A) for pure spherical 
MgSiOa (blue curve) and amorphous carbon (red curve) 
grai ns. Optical c o nstan ts for the silicate dust were taken 
froml.Iager et al.l (2003I1 . and those for the AC dust from 
(IZubko et al.lll99fil b The right panel depicts the wave¬ 
length dependence of the ellipsoidal composite grains. 
The blue curve corresponds to pure MgSiOa grains. The 
black curves give the values of k for the composite grains 
with different amount of AC inclusions (see figure cap¬ 
tion). The green curve corresponds to the mass absorp¬ 
tion coefficient of the composite grains used to fit the 
8515 and 9090 epochs spectra. 

We assumed that the dust radiates at a sin gle tem¬ 
perat ure and used the IDL routine MPFIT dMarkwardtl 
120091 ) to simultaneously fit the SN1987A flux densities for 
the two epochs, with the dust masses and temperatures 
as variables. We experimented with different concentra¬ 
tions of AC inclusions. The results presented here are 
for composite grains with MgSiC >3 and AC volume fill¬ 
ing factors of f s = 0.82 and f c = 0.18, respectively. The 
choice of parameters was driven by the constraints on the 
available mass of Mg and Si in the ejecta. The result¬ 
ing spectra are shown in Figure © and the resulting dust 
masses and temperatures are summarized in Table © 

The mass fraction of MgSiC >3 and AC in the grain are 
given by f s (p s /p) and f c (p c /p), respectively, where p s = 
3.2 g cm 3 and p c = 1.8 g cm 3 are, respectively, the mass 
densities of MgSiC >3 and AC grains, and p = f s p s + f c p c 
is the average density of the grain. If M d is the total 
dust mass, then the mass of Mg, Si, and C locked up in 
the dust is given by 0.24 f s (p s /p) M dl 0.28 f s (p s /p) M d , 
and f c (pc/p) M<j, respectively. For an AC filling factor of 
0.18 we get p = 2.95 g cm 3 , MgSiC >3 and AC mass frac¬ 
tions of 0.89 and 0.11, respectively, and Mg and Si mass 
fractions of 0.21 and 0.25, respectively. The elemental 
masses locked up in the dust are presented in Table [2] 
About half of the MgSiC >3 mass consists of O which, be¬ 
cause of its large abundance in the ejecta (~ 1 M 0 ), does 
not impose any constraints on the mass of dust that can 
form in the ejecta. 

The results show that the late-time spectra require ~ 
0.09 M 0 of Mg, ~ 0.1 M 0 of Si, and ~ 0.046 M 0 of C 
to be locked up in dust. From Table Q] we see that these 
masses are well within their expected yields, alleviating 
the abundance problems of previous models. 

3.2.2. Days 615, and 775 

The most detailed spectra of SN1987A during i ts early 
evolut ionary stages were presented by IWooden et all 
(Il993fl . As mentioned before, the lack of spectral feature 
has been construed as evidence for the lack of silicate 
dust in the ejecta. To test this conclusions we ran mod¬ 
els with pure spherical silicate (MgSiOs) and AC grains, 
in which the mass of Mg and Si in the dust was taken to 
be equal to that derived for epochs 8515 and 9090. For 
this choice of dust masses, or any other sufficiently large 
dust mass, the ejecta is optically thick, and consequently 



A (/im) 



Figure 2. Fits to the observed spectrum with composite grains 
consisting of MgSiOa and AC with volume filling factors of 82 and 
18%, respectively. Detailed description of the fits and tabulated 
masses are given in the text and in Table [2] 

the shape of the dust grains will not affect the blackbody 
spectrum of the ejecta. Our choice of using distinct sili¬ 
cates and AC dust populations is motivated by the dust 
grow th and coagulation models of (ISarangi fc CherchnefB 
I2015H . In these models dust grows primarily by accretion 
at early epochs, after which the abundance of condensible 
elements is sufficiently depleted and coagulation becomes 
the main process for grain growth. 

Thus the only free parameters were the silicate and AC 
dust temperatures. The resulting spectra are shown in 
Figure [3] The red and blue curves are, respectively, the 
carbon and silicate emission components of equation © , 
but only their sum (green curve) is observable. The 
20 /im optical depth of the ejecta is ~ 10 4 and 6000 for 
these two epochs (see also Table©, giving rise to smooth 
featureless spectra. They provide a surprisingly good fit 
to the observations considering the fact that we forced 
each of the silicate and carbon dust to radiate at single 
(but different) temperatures. The figure shows that the 
silicate 9.7 and 18 pm emission features have been in¬ 
ternally absorbed by the silicate-carbon dust mixture, as 
evident from the strong absorption features in the carbon 
dust emission. 

3.2.3. Day II 44 

We selected the 10 and 20 pm photometric data ob¬ 
tained with the Cerro Tololo Inter-American Observatory 
(CTIO) on day 1144 to represent t he observations take n 
during the 935 to 1352 d interval (jSuntzeff et al.l [19911 1. 
For day 1144 we also have a 50 pm upper li mit r e porte d 
by P. M. Harvey, and presented bv iDwek et al.l (|1992f ) . 
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Table 2 

Evolution of Dust Parameters and Ejecta Radius 1 


Epoch (day) 

615 

775 

1144 

8815 

9090 

Pure silicates (MgSiOs) 

M d(AZ©) 

0.40 

0.40 

0.40 



Td(K) 

607 

334 

145 



Pure amorphous carbon 

M d(AZ©) 

0.047 

0.047 

0.047 



T d (K) 

454 

333 

120 



Composite grains z 

M d(AZ©) 




0.42 ±0.07 

0.45 ± 0.08 

T d (K) 




26.3 

25.2 

Elemental abundances in dust 

Mg(M 0 ) 

0.096 

0.096 

0.096 

0.090 

0.095 

Si (M 0 ) 

0.11 

0.11 

0.11 

0.11 

0.11 

C(M 0 ) 

0.047 

0.047 

0.047 

0.045 

0.048 

R e j [10 ib cm] 3 

0.48 

0.61 

0.90 

6.9 

7.1 

L ir (Lq) 4 

3.6 x 10 s 

1.2 x 10 s 

7.0 x 10 3 

310±27 

270±19 

r(A) 5 

9130 

5650 

252 

0.33 

0.33 


1 Entries in bold were held constant. The spectra from the last two epochs were fit with 
ellipsoidal composite grains. Grains in the first three epochs consist of a mixture of spherical 
silicates (MgSiOs), and amorphous carbon (AC) grains. 

2 The composite grains consist of an MgSiOs matrix with AC inclusions occupying 18% of the 
grain’s volume. 

3 The ejecta radius, R e j, at each epoch was calculated for a constant expansion velocity of 
910 km s _1 . 

4 For comparison, on day 1144 the U to 20 ^m (uvoir) luminosity is 3400 L© dSuntzeff et al.l 

emd). 

5 The optical depth, r(A) is given by Equation ([3J), and calculated at around the peak wave¬ 
length of the emission for each epoch, at 20, 20, 50, 200, and 200 fim, respectively. 
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Figure 3. The spectra on days 615 and 775, obtained with a 
combination of pure MgSiOs silicates (blue line) and amorphous 
carbon (AC) grains (red line). Masses of the two components were 
fixed and determined by the masses of Mg, Si, and C locked up 
in dust in epochs 8515 and 9090. The total flux is given by the 
green line. Detailed description of the fits and tabulated masses 
are given in the text and in Table [2] 


X (/xm) 

Figure 4. The spectrum corresponding to the same dust composi¬ 
tions as in epochs 615 and 775, but with temperatures constrained 
to obey the observed upper limit on the IR flux at 50 /xin. Detailed 
description of the fits and tabulated masses are given in the text 
and in Table [2] 

A low resolution (A/AA = 40) 16-30 /xm spectrum of 
the SN was obtai ned with the KAO around day 1151 
dDwek et alJll992 ). Several emission lines were detected, 
and the upper limit set on the continuum intensity is con¬ 
sistent with the CTIO 20 /im detection. Figure 0] presents 
the CTIO data and the 50 /xm upper limit. As for days 
615 and 775, we assumed that the dust giving rise to the 
emission consists of a combination of spherical MgSiOs 
and AC grains with masses identical to this in epochs 615 
and 775, and allowed the temperature to vary. Figure 0] 
shows the fit to the observed spectrum. The dust masses 
and temperatures are presented in Table 0] 

4. CLUMPY EJECTA 

So far we assumed that the ejecta can be approximated 
by a dusty sphere of radius R. However, the early detec- 
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tion of X- and 7 -rays from the SN, and the high ioniza¬ 
tion rate in its envelope, required the macroscopic mixing 
of clumps conta ining radioactive Ni, Co, and F e through¬ 
out the ejecta (|Orahami 11988 : Li et alJll993H . Further¬ 
more, the lack of silicate emission features in the early 
IR spectra (days 615, 775) was attributed to the fact 
that the dust may have formed in optically thick clumps 
dLucv et al.lll989l I1991H . In the following we show that 
the results of our calculations can be simply extended to 
a clumpy medium. 

In the following we will construct a simple model for 
a clumpy ejecta. Since the ejecta is optically thin at 
late epochs (days 8515 and 9090), the model is relevant 
only for the epochs when the homogeneous ejecta was 
optically thick. Using the same dust mass and compo¬ 
sition as in the homogeneous model, we show that the 
clumpy model can reproduce the same IR spectrum at 
each epoch, and is also capable of accounting for the 
UVO region of the SN light curve. 

We will assume the presence of N c identical dusty 
clumps with a dust mass m c = M d /N c , where M d is 
the total dust mass in the ejecta, and radius r c . We also 
assume that the clumps maintain pressure equilibrium 
with their surroundings so that the ratio r c j R, where R 
is the outer radius of the ejecta, remains constant with 
time. We also assume the presence of some dust mass, 
TOj C , in the interclump medium. 

The volume filling factor of the clumps is given by: 

fv = N c x 3 . (7) 


The number of clumps and their radii are constrained 
by the requirement that the projected area of the clumps 
be equal to that of the the homogeneous sphere. This 
requirement reduces to the constraint that N c r 2 = £ R 2 , 
where £ is a factor of order unity that compensates for 
the overlapping of clumps along a given line of sight. 

The optical depth of each clump, r c , is simply related 
to 77, that of the homogeneous sphere, by 


(M d /N c \ _ r h 

4 V 7rr c / K £ 


( 8 ) 


The clumps are therefore optically thick at UVOIR wave¬ 
lengths at all early epochs before day 1144, so that their 
cross section can be taken as equal to the geometrical 
value 7 r r 2 . 

The probability that at photon passes a distance z 
through the ejecta without absorption, and then hits a 
clump in the z and z + dz interval is given by: 


The fraction /a of the projected area of the ejecta that 
is filled by the clumps is then given by the integral: 

f R 

/a = (tuR 2 )- 1 / P(b) 2irbdb (11) 

Jo 

Table [3] summarizes some clump characteristics and 
the fraction of the projected surface of the ejecta cov¬ 
ered by clumps for several values of N c and £. The 
most desirable combination of these parameters is one 
that gives a small volume filling factor for the clumps, 
and a large surface covering factor. The combination 
of {W c , £}={1000, 2} fills both requirements. It give a 
clump filling factor of 0.09, a surface covering fraction 
of 0.82. The value of r c /R for this model is 0.044, con¬ 
sistent with the scales of clumps predicted by numeri¬ 
cal simu lations of the instab ilit ies in core collapse super¬ 
novae (jWongwathanarat et al.ll2015f l. This means that 
in the clumpy model, the integrated brightness of the 
IR emission will be lower by 18% compared to that of 
the homogeneous sphere. This factor is well within the 
uncertainties in the fluxes observed on days 615-1144. 
So all dust masses derived assuming homogeneous ejecta 
are still valid for the clumpy ejecta model. Figure [5] 
presents Monte Carlo simulations of a clumpy ejecta for 
N c = 1000 clumps and a values of £=1.0, 2.0, and 4.0 
(top to bottom). The simulations were performed to 
validate the analytically derived area covering fractions 
listed in Table [3] From top to bottom, the fraction of 
the projected area covered by the clumps is 0.59, 0.82, 
and 0.96, in good agreement with the analytical formula 
presented in Equation (EH). 

The clumpy model must allow for the escape of UVO 
photons, and any selective extinction in the ejecta. In the 
followin g; we use th e mega -grain approximation described 
bv lYarosi fc I) weld (1 999h to calculate the escape of these 
photons from the clumpy medium. Equation (59) in their 
paper can be simplified by assuming that the volume 
filling factor, /y, of the clumps is small, that most of 
the dust resides in the clumps, and that the clumps are 
solid spheres with geometrical cross section nr 2 . The ra¬ 
dial optical depth, tr( A), of the ejecta is then given by 
the sum of two terms. The first is wavelength indepen¬ 
dent and caused by the absorption in the optically thick 
clumps, and the second is the wavelength dependent ab¬ 
sorption by dust residing in the interclump medium. 

tr{\) = n c nr 2 R + pi C k(A) R (12) 

= (3/4) [£ + ( mi C /nR 2 ) k( A)] 


P(z)dz = exp(-z/£) — (9) 

where t = (n c 7rr 2 ) -1 is the photon mean free path 
through the ejecta, and n c is the number density of 
clumps. A photon incident on the sphere within an an¬ 
nulus 2irbdb, where b is the projected distance from the 
center of the sphere, will strike a clump with a probabil¬ 
ity given by: 

f 2Ry/l-(b/Ry 1 

mdb = J o exp (-z/£)- ( 10 ) 

= 1 -exp (-2Ry/1 - ( 6 /i?) 2 ) 


where pi C and rrii C are, respectively, the density and mass 
of the dust in the interclump region, and V is the volume 
of the ejecta. 

At early times, the IR emission accounts for half the 
energy deposited in the ejecta by the radioactive decay 
of 56 Co. So the unabsorbed UVO luminosity can be at 
most equal to that of the IR. The observed UVO flux 
constitutes about 5-10% of the IR emission fr om the SN 
(jLucv et all 119891 11991b iWesson et al.1 120151 Figures 1 
and 2). So at least 5 to 10% of the unabsorbed UVO 
emission must escape the interclump medium. These, 
and observations of neb ular emission lines from the ejecta 
(jjerkstrand et al.ll 2011 lf . constrain the mass of dust that 
can be present in this medium. 
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Table 3 

Surface covering factor 
for several clump models 


N c 

e 

fv 

r c /R 

/a 

100 

1 

0.10 

0.10 

0.60 


2 

0.28 

0.14 

0.82 


4 

0.80 

0.20 

0.94 

1000 

1 

0.03 

0.03 

0.60 


2 

0.09 

0.044 

0.82 


4 

0.25 

0.063 

0.94 


Figure 5. Monte Carlo simulations of clumpy ejecta characterized 
by N c = 1000 and $=1.0, 2.0, and 4.0 (top to bottom). The column 
on the left depicts the dust column density. The brighter regions 
depict regions of overlapping clumps. The right column depicts 
the projected surface area of all the clumps. The surface covering 
factors, f a for this simulation are listed in Table 3. 

For radial optical depths > 1, the escape probabil¬ 
ity from the ejecta is approximately given by 3/4r. 
For £ = 2, the first term in Equation (1121) contributes 
3£/4 = 1.5 to t Ri giving P esc « 0.5. For P esc to be larger 
than 0.05, the optical depth must be smaller than ~ 15, 
setting a limit of r < 13 on the second term in eg. (fl^l) . 
The mass absorption coefficient of the dust at 0.55 fim is 
about 10 4 cm 2 g _1 for a mixture of silicate-carbon grains 
(Figure 1), so the second term becomes 

(3/4) ( m ic /nR 2 ) k(V) « 2 x 10 5 m ic (M 0 ) < 13 (13) 

where we assumed an ejecta radius of 0.5 x 10 16 cm. So 
the interclump dust mass is limited to rrii C < 10 -4 Mq, 
allowing UVO emission that is not absorbed in the 
clumps to escape the ejecta. 

In summary, the clumpy model can reproduce most 
of the early IR emission from the ejecta with the same 
mass of dust as in the homogeneous model. The clumps 
are optically thin at late times, so that the mass derived 
for days 8515 and 9090 is unaltered as well. The mass 
of dust in the interclump region is small compared to 
that in the clumps, allowing for the escape of nebular 
and continuum UVO emission from the ejecta. Because 
of the small mass of dust in the interclump medium, it 
makes no significant contribution to the IR emission from 
the clumps. 


Table 4 

Inferred 44 Ti Yield from UVOIR and FIR Observations 
of SN1987A 


M( 44 Ti) 
(10- 4 Mq) 

Observations 

Epoch 

(days) 

Reference 1 

3.1 ±0.8 

hard X-rays 

3600-4000 

[i] 

1.5 

UVOIR lines 

~3000 

[2] 

0.55 ±0.17 

UVOIR 

<1900 

[3] 


extrapolated V 

<4300 


1.6 

far-IR dust 

8500-9100 

this work 


1 References: [1 Grcbcncv et al. (2012); [2] 

IJerkstrand et a.Tl 1120 M l: [3] Seitcnza hl et al.l 120141) 


5. THE EVOLUTION OF THE IR LUMINOSITY AND THE 
44 TI ABUNDANCE IN THE EJECTA 

5.1. The Evolution of the IR Luminosity 

The IR luminosity at each epoch is simply given by 
the integral of the spectrum, assuming a source distance 
of of 50 kpc. The results are given in Table [2] At early 
times, the IR luminosity follows the 56 Co decay curve, 
but is consistent with being constant on days 8515 and 
9090. 

The energy absorbed and re-radiated by the dust could 
be either generated by sources external to the ejecta, or 
by internal energy sources. External energy sources in¬ 
clude the X- ray and UVO emis sion generated by the re¬ 
verse shock dFrancc et al.ll2015l) . and/or the diffuse inter¬ 
stellar radiation field and nearby stellar sources. Internal 
energy sources include a potential embedded pulsar or 
decaying radioactive isotopes. 

5.2. External Energy Sources 
5.2.1. The reverse shock 

Prediction that the reverse shock generated by the 
interaction of the SN blast wave with the circumstel- 
lar ring would be manifested in increased line emission 
from the ejecta were confirmed by UVO imaging spec¬ 
troscop y of SN1987A with the Hubble Space Telescope 
( HST) (IFrance et al.ll201ll and references therein). Dur¬ 
ing the period of interest here (~ day 9700), the X-ray 
- UV emission of the ER and the reverse shock over 
the 0.0 1 -8 keV en ergy range is ~ 1460 L & (day 9885, 
IFrance et al.l 1201 5:, apparently sufficient to power the 
IR emission. However the solid angle subtended by the 
ejecta (~ 2 x 10 17 x 1 x 10 17 cm), as seen from the equa¬ 
torial ring (ER) (radius = 6.1 x 10 17 cm) is only 0.158 sr. 
So even if all this incident radiation traverses the dust- 
free external H-layers and gets absorbed by the dust, only 
~ 1.3% of this luminosity can be absorbed and reradiated 
by the ejecta. 





























5.2.2. Diffuse and discrete stellar sources 

Other external heating of the ejecta can be supplied 
by the interstellar radiation field (IRSF). However, if the 
energy density is the same as in the solar neighborhood 
(jMathis et al.1 1 19831 iGalliano et al.l 120111) , then the in¬ 
cident flux is 55 times weaker than the (unattenuated) 
X-ray flux. 

Nearby Stars 2 and 3 are B2 I II stars with lumi nosi ties 
of ~ 2 x 10 4 and ~ 1 x 10 4 L 0 (|Walborn et al.lll993l) . If 
these stars are at the same line of sight distance as SN 
1987A, then the stars can provide fluxes as large as 1 and 
1.5 times (respectively) the flux of the X-rays. 

Altogether, the X-ray emission can at most provide 
1460 x 0.158/47T = 18 L 0 . The local ISRF can provide 
only 0.3 Lq, and stars 2 and 3, at most 18 and 27 Lq, 
respectively. So at most, ~ 63 L 0 of the total IR emission 
of ~ 300 L 0 , can be supplied by external energy sources. 

5.3. Internal Energy Sources 
5.3.1. Compact central object 

One internal source of energy would be from accretion 
onto or spin down of a pulsar. Searches for a compact 
remnant have been unsuccessful (except possibly Mid- 
dleditch et al. 2000), and hav e set low lumi nosi ty limits 
on a compact source (< 2L 0 ) dGraves et al.ll2005l) . How¬ 
ever, this could be higher if the extinction is far larger 
than was assumed. 


5.3.2. Decaying radioactive elements 

During the first two years after the explosion, the SN 
light curve was powered by the energy released in the 
radioactive decay chain 56 Ni—> 56 Co—» 56 Fe. Models 
show that ~ 0.07 M 0 of 56 Ni were required to fit the 
early U VOIR light curve, b e fore dust has formed i n the 
ejecta (TWooslev et al.1119881 ISeitenzahl et al.l 12014 and 
references therein). After about 2000 days the SN1987A 
light curve is predominantly powered by the decay of 
44 Ti (ji /2 = 58.9 yr). Efforts to determine the amount 
of 44 Ti required to power the late time light curve have 
focused on hard X-rays, nebular lines, and UVOIR ob- 
servations of the SN. A detailed summary is given by 
ISeitenzahl et al.l (|20l4 , with highlights given in Table 0] 
of this paper. Until day ~ 1800 the luminosity was deter¬ 
mined from the UVOIR observations, which extend out 
to the 20 /Mil band. However, most of the dust emission 
after day ~ 1100 is emitted at longer wavelengths (see 
Figure 01 so that he UVOIR luminosity may be missing 
a significant fraction of the total energy emitted by the 
ejecta. 

Figure [6] compares model light curves to the IR lu¬ 
minosities (yellow diamonds). The 56 Ni, 57 Ni, 55 Co, 
and 60 Co, and masses were taken to be equal to 0.071, 
4.1 x 10~ 3 , 9.2 x 10~j and 4.5 x 10~ 8 M 0 , respec¬ 
tively (ISeitenzahl et al.112014 . The figure shows that 
at early epochs (days 615, 775), the IR emission ac¬ 
counts for about half the radioactive energy deposited 
by the decay of 56 Co in the ejecta. At later epochs, the 
IR luminosity can provide an important measure of the 
amount of 44 Ti created in the explosion. The different 
curves correspond to different masses of 44 Ti with val¬ 
ues of 0.55 x 10 -4 ( Seitenzahl et al.l 12014 red curve), 
and 3.1 x 10~ 4 dGrebenev et al.ll2012t blue curve). Also 



Figure 6. The energy released by the decay of radioactive ele¬ 
ments in the ejecta as a function of time. All luminosities were 
calculated for a 56 Ni abundance of 0.071 Mq, but with different 
44 Ti abundances, spanning the range of values in Tabic HI The yel¬ 
low diamonds represent the IR luminosities emitted by the dust at 
the different epochs (Table |2j. The black curve represents the best 
fit to the luminosities on days 8515 and 9090. The open diamond 
represents the T-ReCS observations on dav 6067 IBouchet et al .1 
12004 . 

shown in the figure is the luminosity of 390 L 0 on day 
6067, derived from the T-R eCS observations of the SN 
ejecta (IBouchet et al.ll2004 open diamond). The black 
curve provides the best fit to the IR luminosities derived 
in this paper and requires M( 44 Ti)= 1.6 x 10~ 4 Mq. 

This val ue is in good agreement w ith the 1.5 x 10 -4 M 0 
derived by l.lerkstrand et al.l 1201 If) from the analysis of 
UVO emission lines in the ejecta around day 3000. How¬ 
ever, in their standard model they adopted a dust optical 
depth of only 1. The IR optical depth on that day could 
be around 10 at the peak of the IR emission, and sig¬ 
nificantly higher at UVO wavelengths. The 44 Ti mass 
inferred from the energy required to power the nebular 
emission lines should therefore be increased to account 
for the effect of dust absorption. 

On the other hand, the 44 Ti mass inferred from the IR 
emission alone should also be modified: increased to ac¬ 
count for any UVO emission, and decreased to account 
for the contribution of external heating sources to the 
emission. The mass should be definitely increased by 
~ 2%, to account for the UVO luminos ity of ~ 7 L 0 ob¬ 
served on day 8328 (lLarsson et al.ll2013l Figure 5), and at 
most decreased by ~ 20% to account for the possible con¬ 
tribution from external heating sources. So the IR obser¬ 
vations provide a robust lower limit of ~ 1.3 x 10 -4 M 0 
on the mass of 44 Ti synthesized in the explosion. 

6. SUMMARY AND ASTROPHYSICAL IMPLICATIONS 

In this paper we presented a new scenario for the evolu¬ 
tion of dust in SN1987A. In this scenario, the observed IR 
spectrum on days 8515 and 9090 arises from ~ 0.45 M 0 
of composite ellipsoidal grains, consisting of a mixture 
of MgSiCU and amorphous carbon with volume filling 
factors of 82 and 18%, respectively. I n contrast with 
previous models dMatsuura et al.l 120151 ). this model re¬ 
quires only ~ 0.1 M 0 of Mg and a similar mass of Si to 
be locked up in dust, consistent with the post-explosive 
mass of these elements in the SN ejecta. The mass of 
carbon locked up in du st is only 0.04 Mq, in c ontrast to 
the 0.5 M & of C in the lMatsuura et all (120151 ) models. 

We also show that the same mass of Mg, Si, and C 
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could have already been locked up in the form of spherical 
MgSiC >3 and AC grains since day 615 in IR-optically thick 
clumps. The low mass of dust and the upper limit on the 
silicate dust in previous models were inferred from the 
absence of the 9.7 and 18 ^m silicate emission features, 
assuming optically thin clumps. However, we show that 
the absence of the silicate absorption features at the early 
epochs is the result of the large optical depth and self¬ 
absorption in the clumps. Figures[2j[3l and|4]and Table[2] 
show, respectively, the fit to the observed spectra, and 
the dust parameters corresponding to the fits. 

The model presented in this paper t hus s ol ves t hree 
problems that plague pr evious models dMatsuura et al.1 
120151 iWesson et al.ll2015l and refernces therein): 

1 . the excessive mass of silicate or AC dust required 
to fit the late-time spectrum 

2 . the need for cold accretion to grow the grains from 
~ 10" 3 to ~ 0.5 Mq 

3. the absence of SN-condensed silicate dust 

The results presented in this paper a re si g nificantly 
differen t fro m the conclusi o n rea ched by lErcolano et al.1 
(|2007af) and IWesson et al.1 (|2015T) , who firmly ruled out 
the possible presence of more than ~ 3 x 10 -3 Mq of dust 
at the early epochs before ~ day 1200. The main differ¬ 
ence between our respective models is that their mod¬ 
els were constructed to simultaneously fit the observed 
UVOIR as well as the IR dust emission. Without this 
constraint, our models can fit the mid- to far-IR emis¬ 
sion at early times with dust having higher temperatures 
and densities, being optically thick through far-IR wave¬ 
lengths. At the high densities prevailing at early epochs, 
the dust is thermally coupled to the gas, and more likely 
to be collisionally heated by the ambient hot gas, rather 
than radiatively heated by UVO photons. In either case, 
the ultimate heating sources are the decaying radioactivi¬ 
ties in the ejecta. This scenario seems to be confirmed by 
the fact that the late time IR spectra, when the ejecta 
is optically thin, were fitted with a very narrow range 
in dust temperatures. Such temperature distributions 
can be more readily explained if the dust is collision- 
ally heated in the ejecta. In our scenario, and UVOIR 
emission solely arises from the sources in the interclump 
region, where the optical depth is very low. Only ~ 10% 
of the energy is emitted at these wavelengths after the 
dust formation epoch. 

The IR observations provide a robust lower limit of 
~ 1.3 x 10~ 4 Mq on the mass of radioactive 44 Ti that 
has formed in the explosion, with a more likely value of 
1.6 x 10 -4 Mq, depending on the contribution of external 
sources in heating the dust. 

The proposed scenario requires the very rapid evolu¬ 
tion of dust in the SN ejecta. Theoretical calculations 
show that only ~ 3 x 10~ 3 of dust can form by day 
700 assuming a smooth ejecta, a value that increases to 
~ 0.08 Mq for clumpy ejecta (e.g. iSarangi fe Cherchnefll 
120151! . Further increase in the dust mass can be achieved 
by increasing the density in the clumps, and by the inclu¬ 
sion of molecular and dust-continuum cooling, processes 
that were ignored during the dust formation epoch. 

Finally, our model cast doubts on interpretations for 
the existence of an evolutionary trend in the dust mass 
from young SNe to older remnants. Infrared observations 
of young SNe need to carefully determine the origin of 


the emission, excluding any possible contributions from 
IR echoes and shocked interstellar dust. An extensive 
review of the observations of the IR emission from SNe, 
and their rol e as dust s o urces i n the early universe was 
pr esented bylGall et all (|2011bH . In p articular. F i gure 2 
of lGall et all (|2011bh and Figure 4 in I Gall et al.1 (|2014H 
show that the inferred dust mass in SNe that are less than 
~ 3 yr old, is between ~ (0.5 — 50) x 10 -4 Mq, whereas 
the mass of SN condensed dust inferred from observa¬ 
tions of young unmixed remnants, such as Cas A and 
the Crab Nebula i s about 0.1 M ^ (lBarlow_et_aL| EMS: 
Arendt et al.1 12014 I Gomez et all 120121 : Temim fc Dwekl 
2013f) . This “trend” could be interpreted as an evo¬ 
lutionary one, suggesting the growth of dus t grains i n 
the eje cta with time. However, as noted by iGall et ahl 
(!2011bjh the low dust mass in young SNe could also be 
an observational selection effect, since it is inferred from 
only near-IR observations, leaving a large part of the dust 
spectrum unexplored. In this paper we showed that even 
when the rnid-IR region of the spectrum is observed, op¬ 
tical depth effects can result in a large underestimate of 
the dust mass in the ejecta. 

We thank Arkaprabha Sarangi for helpful discussions, 
and the referee for useful criticism that led to the im¬ 
provement of the manuscript. This work was supported 
by NASA Astrophysical Data Analysis Program 13- 
ADAP13-0094. 
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